Light Quark Masses 99 * 



Stephan Narison a 

a Laboratoire de Physique Mathematique, Universite de Montpellier II Place Eugene Bataillon, 34095 - 
Montpellier Cedex 05, France and KEK Tsukuba, Ibaraki, 350801 - Japan. 
E-mail: narison@lpm.univ-montp2.fr 

I give a short historical and a critical review of the determinations of light quark masses from QCD at dawn of the 
next millennium. QCD spectral sum rules combined with ChPT give, to order af, the world average for the running 
masses: m s (2 GeV) = (118.9 ± 12.2) MeV, m d (2 GeV) = (6.3 ±0.8) MeV, m„(2 GeV) = (3.5 ±0.4) MeV and the 
corresponding values of the invariant masses given in Eq. 24. Lower and upper bounds derived from the positivity 
of spectral moments are presented in Tables 2 and 3. For a comparison, we critically review the recent lattice results 
(section 8 and Table 5) and attempt to deduce the present QCD grand average determination (Table 6): m s (2 
GeV)= (110.9±8.8) MeV, to be used with a great care. Then, we deduce the value: B^ 2 - 0.45(resp. 0.32)Bg /2 ~ 
1.6±0.4 (resp. 1.1±0.3) and the lower bound 1.1±0.2 (resp. 0.7±0.1), for the combination of the penguin operators, 
governing the CP-violating parameters e' /e without (resp. with) the inclusion of the final state interaction effects. 
The result signals a possible deviation from the leading 1/N prediction by about 1 ~ 3<r, which should be tested 
using accurate non-perturbative calculations. 



1. INTRODUCTION 

One of the most important parameters of the 
standard model and chiral symmetry is the light 
quark masses. Indeed, they are useful for a much 
better understanding of the realizations of chiral 
symmetry breaking and for some eventual 

explanation of the origin of quark masses in uni- 
fied models of interactions . Within some pop- 
ular parametrizations of the hadronic matrix ele- 
ments H , the strange quark mass can also largely 
influence the Standard Model prediction of the 
CP violating parameters e'/e, which have been 
mesured recently ||. However, contrary to the 
leptons, where the physical masses can be identi- 
fied with the pole of the propagator the quark 
masses are difficult to define because of confine- 
ment. Instead, they can be treated as coupling 
constants of the QCD Lagrangian, where the no- 
tion of the running and invariant masses, which 



*Review talk given at the QCD 99 Euroconference (Mont- 
pellier 7-13th July 1999) and plenary talk given at the 
QCD Confinement 2000 (Osaka 7-10th March 2000). 
2 For a first explicit definition of thenerturbative quark 
pole mass in the MS-scheme, see |7J (renormalization- 
scheme invariance) and ji| (regularization-scheme invari- 
ance) . 



are renormalization scheme and scale dependents, 
has been introduced In practice, these masses 
are conveniently defined within the standard MS- 
scheme. In addition to the determination of the 
ratios of light quark masses (which are scale inde- 
pendent) from current algebra jjj, and from chiral 
perturbation theory (ChPT), its modern version 
||, a lot of effort reflected in the literature JTo[ 
has been put into extracting directly from the 
data the running quark masses using the SVZ 
QCD spectral sum rules (QSSR) LEP 
experiments [Q and lattice data |l3]. In this 
talk, I shall review the different determinations 
from these QCD approaches, by emphasizing the 
historical developments of the field. 

2. RUNNING AND INVARIANT LIGHT 
QUARK MASSES IN QCD 

It is convenient to introduce the dimcnsionlcss 
coupling Xi(i>) = mi(is)/is, where v is the renor- 
malization scheme subtraction constant. The 
running quark mass is a solution of the differ- 
ential equation: 

dx ' 

— l - = (l+7(a,))s i (t) : Xi(t = 0) = Xi (y) . (1) 



In the AfS'-scheme, its solution to order (a s 
a s /n) is: 

ffh(y) = m i (-p 1 a s (v))-"> l/f31 |l 



+ 



02 ( 71 72 \ 



1 


[1 


( 7i 


\ 2 

72 \ 


Pi 


( 7i 


72 


2 




Ui 


02 J 


Pi 


\Pi 


02 



03 f 7i _ 73 

/3l Ul 08 



1.95168^ ^ , 



(2) 




where the term comes from ]S ; 7^ are the 
C(a*) coefficients of the quark-mass anomalous 
dimension, which read for three flavours: 



Figure 1. m s /m d versus m u /md from |Q . 



7i 



72 = 91/12, 73 = 24.8404 



(3) 



The invariant mass rhi has been introduced for 
the first time by || in connection with the analy- 
sis of the breaking of the Weinberg sum rules by 
the quark mass terms in QCD. 

3. RATIOS OF LIGHT QUARK MASSES 

The ratios of light quark masses are well- 
determined from current algebra pL and ChPT 
0. In this approach, the meson masses are ex- 
pressed using a systematic expansion in terms of 
the light quark masses: 



M 2 K+ 



(m u +m d )B + 0(m 2 ) + 
(m u +m s )B + 0(m 2 ) + . 
(m d + m s )B + 0(m 2 ) + . 



(4) 



where B = — (ipijj) / ffc from the Gell-Mann, 
Oakes, Renner relation fig] : 



2 f 2 



-(m u + m d )(ipijj) + 0(m 2 



(5) 



However, only the ratio, which is scale indepen- 
dent can be well determined. To leading order in 



m llTjfr 

m u ^ M 2 + - 

m d ~ M 2 + 4 

m. ^ -M* + 

m d ~ M 2 - 



M 2 K0 - M 2 K+ 



M 2 K0 



M| 



0.66 



20 



(6) 



Including the next order + electromagnetic cor- 
rections, the ratios of masses are constrained on 
the ellipse: 



m d 



1 / m s 
Q 2 \m d 



(7) 



where: Q 2 ~ (m 2 s - m 2 )/{m 2 d - ml) = 22.7 ± 0.. 



using the value of the 77 



from the PDG 



average [|10| , though this value can well be in the 
range 22-26, to be compared with the Dashen's 
formula @ of 24.2; rh = (\/2)(m u + m d ). In Fig. 
1, one shows the range spanned by R = (m s — 
rh)/(m d — m u ) and the corrections to the GMO 
mass formula A M : Mf = (1/3)(4M|- - m 2 r )(l + 
Am)- The Weinberg mass ratio |Q is also shown 

3 In Generalized ChPT, the contribution of the m 2 -term 
can be as large as the m one Jli| , which modifies drasti- 
cally these ratios. 



which corresponds to the Dashen's formula and 
R ~ 43. At the intersection of different ranges, 
one deduces [|l7f : 



m„ m, 
— = 0.553 ±0.043 , — 

m d m d 

- f_ =24.4 ±1.5. 

(m d + m u ) 



18.9 ±0.8, 



(8) 



The possibility to have a m u = advocated in 
20 appears to be unlikely as it implies too strong 
flavour symmetry breaking and is not supported 
by the QSSR results from 2-point correlators of 
the divergences of the axial and vector currents, 
as will be shown in the next sections. 

4. QCD SPECTRAL SUM RULES 

4.1. Description of the method 

Since its discovery in 79 jll), QSSR has proved 
to be a powerful method for understanding the 
hadronic properties in terms of the fundamental 
QCD parameters such as the QCD coupling a s , 
the (running) quark masses and the quark and/or 
gluon QCD vacuum condensates. The description 
of the method has been often discussed in the lit- 
erature, where a pedagogical introduction can be, 
for instance, found in the book jl2). In practice 
(like also the lattice) , one starts the analysis from 
the two-point correlator: 

i> H {q 2 ) = i [ d A x e ic > x (0\TJ H (x) (JH(0)) f |0) ,(9) 



built from the hadronic local currents Jh{x), 
which select some specific quantum numbers. 
However, unlike the lattice which evaluates the 
correlator in the Minkowski space-time, one ex- 
ploits, in the sum rule approaches, the analyt- 
icity property of the correlator which obeys the 
well-known Kallen-Lehmann dispersion relation: 



<M? 2 ) 



dt 



t-q 2 



- Im^(i) + ..., (10) 



where ... represent subtraction points, which are 
polynomials in the g 2 -variable. In this way, the 
sum rule expresses in a clear way the duality 
between the integral involving the spectral func- 
tion ImipH (t) (which can be measured experimen- 
tally), and the full correlator tpniq 2 )- The latter 



can be calculated directly in the QCD Euclidean 
space-time using perturbation theory (provided 
that — q 2 + m 2 (m being the quark mass) is 
much greater than A 2 ), and the Wilson expan- 
sion in terms of the increasing dimensions of the 
quark and/or gluon condensates which simulate 
the non-perturbative effects of QCD. 

4.2. Beyond the usual SVZ expansion 

Using the Operator Product Expansion (OPE) 
PH], the two-point correlator reads: 



1 



D=0,2,. 



(q 2 ) 



2\D/2 



C(q 2 ,v)(0(v)), 



dimO—D 



where v is an arbitrary scale that separates 
the long- and short-distance dynamics; C are 
the Wilson coefficients calculable in perturbative 
QCD by means of Feynman diagrams techniques; 
{0(v)) are the quark and/or gluon condensates 
of dimension D. In the massless quark limit, one 
may expect the absence of the terms of dimension 
2 due to gauge invariance. However, it has been 
emphasized recently pH that the resummation of 
the large order terms of the perturbative series, 
and the effects of the higher dimension conden- 
sates due e.g. to instantons, can be mimiced by 
the effect of a tachyonic gluon mass A which gen- 
erates an extra D = 2 term not present in the 
original OPE. Its presence might be understood 
from the analogy with the short distance linear 
part of the QCD potential []. The strength of this 
short distance mass has been estimated from the 
e+e- data to be lElES: 



—A 2 



-(0.06 - 0.07) GeV 2 



(11) 



which leads to the value of the square of the 
(short distance) string tension: a ~ — |a s A 2 ~ 
[(400 ±20) MeV] 2 in an (unexpected) good agree- 
ment with the lattice result |25|] of about [(440 ± 
38) MeV] 2 . In addition to Eq. 5, the strengths 
of the vacuum condensates having dimensions 
D < 6 are also under good control, namely: 

• (ss)/(dd) ~ 0.7 ± 0.2 from the meson pj 
and baryon systems Pq[; 



4 Some evidence of this term is found from the lattice anal- 
ysis of the static quark potential [ p2| , though the extrac- 
tion of the continuum result needs to be clarified. 



(a s G 2 ) ~ (0.07 ± 0.01) GeV 4 from sum 
rules of e + e~ — > 1=1 hadrons p^] and heavy 
quarkonia [p7|-p9| , and from the lattice (30) ; 

g(i>\ a /2a^G a ^) ~ (0.8±0.1) GeV 2 W), 
from the baryons pi 26 1 and the heavy-light 



mesons 32 
\2 , 



5.8 x 10" 4 GeV b from e+e" 



/ = 1 hadrons 23 1; 



• g 3 (G 3 ) w 1.2 GeV 2 (a s G 2 ) from dilute gaz 
instantons (3^] . 

4.3. Spectral function 

In the absence of the complete data, the spectral 
function is often parametrized using the "naive" 
duality ansatz: 



- 1mip H (t) =s 2Mj?f 2 H S{t - M% 

IT 

"QCD continuum" x fl(t - t c ) , 



(12) 



which has been tested |12| using e + e~ and r- 
decay data, to give a good description of the spec- 
tral integral in the sum rule analysis; fn (ana- 
logue to /„■) is the hadron's coupling to the cur- 
rent ; In is the dimension of the correlator; while 
t c is the QCD continuum's threshold. 

4.4. Form of the sum rules and optimiza- 
tion procedure 

Among the different sum rules discussed in the 
literature we shall be concerned with the: 
• Laplace sum rule (LSR) §I||||28j: 



Cn(r) = 



1 

dt t n exp(-tr) - Imip H (t) ■ (13) 

7T 



The advantage of the Laplace sum rules with re- 
spect to the previous dispersion relation is the 
presence of the exponential weight factor, which 
enhances the contribution of the lowest resonance 
and low-energy region accessible experimentally. 
For the QCD side, this procedure has eliminated 
the ambiguity carried by subtraction constants, 
arbitrary polynomial in q 2 , and has improved the 
convergence of the OPE by the presence of the 
factorial dumping factor for each condensates of 
given dimensions. As one can notice, there are "a 
priori" two free external parameters (r, t c ) in the 



analysis. The optimized result will be (in princi- 
ple) insensitive to their variations. In some cases, 
the i c -stability is not reached due to the too naive 
parametrization of the spectral function. One can 
either fix the t c - values by the help of FESR (lo- 
cal duality) or improve the parametrization of the 
spectral function by introducing threshold effects 
fixed by chiral perturbation theory, in order to 
restore the ^-stability of the results. The results 
discussed below satisfy these stability criteria. 
• r-like sum rules 131- 



Mi 



dt f 



1 



t 

m7 



IrmMt) ,(14) 



The advantage of the r-like sum rule is the pres- 
ence of the threshold factor which gives a zero 
near the real axis where QCD is not expected to 
be applicable. Optimal results should be insen- 
sitive to the changes of M T) and to the values of 
the degrees (m, n) of the moments. 

• Finite Energy Sum Rule (FESR) |§-[|l) : 

M„(t c ) = f dtt n - Jmtp H (t) , (15) 
Jo t 

The advantage of the FESR is the separation (to 
leading order in a s ) of the terms of given di- 
mensions, which gives a set of local duality con- 
straints. However, unlike the two formers, FESR 
is sensitive to the high-energy tails of the spec- 
tral integral and needs an accurate treatment of 
this region, in order that the optimal results are 
insensitive to the changes of t c . 

5. UP AND DOWN RUNNING MASSES 

5.1. Pseudoscalar sum rules 

• Values of (m u + rfld) have been extracted for 
the first time in jl2|,[34j using the sum rule of the 
2-point correlator associated to the pseudoscalar 
current: 



9 M A M (x) = (mj + rrtj) : u(ij 5 )d 



(16) 



The analysis has been improved (or disproved) 
later on by many groups by the 

inclusion of higher order terms or/ and by a more 
involved parametrization of the spectral function 
(threshold effects, ChPT,...). However, this chan- 
nel is quite peculiar due to the Goldstone nature 



Table 1 

(m„ +rnrf)(l) in MeV updated to order 



(m u + m d )(l) 



Sources 



Authors 



Table 2 

Lower bounds on m Uj( j,s(2) in MeV to order a s 



Observables 



Sources 



Authors 



Best estimate 

12.8 ± 2.5 

12.1 ±2.4 
12.6 ±3.2 



Others 
13.2 ±4.4 

15.6 ±3.4 



ir +ChPT BPR95 |6|]4 

BPR+ CNZ99 H 

tachyonic gluon 

(^V) DN98* H 

N, B*-B LSR 
D -> iT*Zz/ LSR 



7r±moments 

LSR 

tt + tt' NWA 



Y97 ( 
SN89 



rn u + m d 
8 

7.3 
7. 



7T, (7 LRT97JM 
7T Y97[|0f 
(^V)±GMOR DN98T53] 



m d - m u 
1.5 



Y97N50 



100 
104 
90 



K 
K 

(V^>)+ChPT 



LRT97p[ 
Y97pf 
DN98l5ll 



13.1 ± 1.5(stat) ± 1.3(syst) Average 



* to order a s and not included in the average. 



of the pion, where the value of the sum rule scale 
(1/r for Laplace and t c for FESR) is relatively 
large of about 2 GeV 2 compared with the pion 
mass, where the duality between QCD and the 
pion is lost. This implies an important role of 
the higher states (radial excitations or/and the- 
oretical parametrizations of the spectral function 
above the 3-7T threshold) in the analysis, and then 
led to some controversial results, which hopefully 
can be cured by the presence of the new 1/q 2 
|2l|]24f | due to the tachyonic gluon mass, which 
enlarges the duality region to lower scale and 
then minimizes the role of the higher states into 
the sum rule. The errors due to the QCD part 
of the sum rules, which is now known to order 
a^, are much less than from the parametrization 
of the spectral function. Among the available 
results, we consider that the best estimates of 
(m u +fr%d) from this channel come from jl6| (ir+ 
ChPT parametrization of the 3n continuum) and 



from |24| (inclusion of the tachyonic gluon mass 
into the analysis of Q). The result of |46|| to 
order c? s has been extended to order by [47|. 
Also the result of [|6| updates the one in [45|. 
Ref. |Q uses the positivity of the higher state 
contributions plus the moment inequalities to or- 
der a s . In Q], one treats the it' in a Narrow 
Width Approximation (NWA), while the QCD 
expression is to order a 2 s . These different deter- 
minations are quoted in Table 1, after including 
into these published results the perturbative con- 
tributions of the known a^-term The effect 
of this term is quite small as the PT series con- 
verges quite well at the sum rule working region of 
about 1.5 GeV. Indeed, the PT expression of e.g. 
the Laplace sum rule normalized to (m u +md) 2 
reads: 

C - 1 + 4.82a s + 21.98a 2 + 53.14^ + C«).(17) 

One can notice the good consistency of the re- 
sults from the different forms of the sum rules in 
the pion channel. We give in Table 1 the aver- 
age of these updated determinations, where we 
have added an extra 10% error which takes into 



Table 3 

Upper bounds on m U! d,s(2) in MeV from |52 
(resp.[^5)) to order a s (resp. a 3 ) 



Observables 



Sources 



Authors 



m u + m d 
11.5 

m s 
148 
147 ±21 



(^V)+GMOR DN98p| 



(^)+ChPT DN98Q 
: + e-+r-decay SN99f35[ 



account the systematics of the approach^]. 

• Lower bounds for (m u + mj) based on mo- 
ments inequalities and the positivity of the spec- 
tral functions have been obtained, for the first 
time, in (4^J|^]. These bounds have been red- 
erived recently in [|^,[5(| to order a s []. Their op- 
timal values quoted in Table 2 exclude the low 
value of about 6 MeV given by Jjsfl . 

5.2. Scalar sum rules 

• Lower bounds on (ma — fn u ) have been ex- 
tracted for the first time in |54| using the sum rule 
of the 2-point correlator associated to the scalar 
current: 



d, j y tl {x) = (m d - m u ) : d(i)i 



(18) 



which is sensitive to leading order to the quark- 
mass difference. The analysis has been extended 
later on by many authors |jl^,^5|,|o| . However, 
the analysis relies heavily on the less controlled 
nature of the a (980), where its qq nature appears 
to be favoured by the present data Q . In the / = 
channel, the situation of the it-it continuum is 
much more involved due to the possible gluonium 
nature of the low mass and wide a meson 1 56 57], 



5 The estimation of the systematic error is based on the 
(un)ability of the method for reproducing the hadron 
masses and couplings |llll3|. 

6 The inclusion of the of^ term will decrease by about 10% 
the strength of these bounds, which is within the expected 
accuracy of the result. 



|55| . Instead, one can also use these sum rules 
the other way around, i.e. by using the values 
of the quark masses from the pseudoscalar sum 
rules and their ratio from ChPT, in order to test 
the nature of these resonances Jl2],0 . 

5.3. Direct extraction of {4>ij}) 
• The chiral condensate can be directly ex- 
tracted from the sum rules (nucleon, B*-B split- 
ting, vector form factor of D* — > K*lv), which are 
particularly sensitive to it and to the mixed con- 
densate {4>a^{X a /2)Gl v il)) = M 2 (V^). A global 
fit from these different channels gives, to order 
a s , the running condensate value at 1 GeV 152] : 



0.6 < (V^)/[-225 MeV] 3 < 1.5, 



(19) 



a result also recovered by the lattice |33[ . 
• Lower and upper bounds on the light quark 
masses given in Tables 2 and 3, can be obtained 
by transforming this result using the PCAC re- 
lation in Eq. 5, and the positivity of the 0(m 2 ) 
term. These results are independent on how chi- 
ral symmetry is realized (ChPT or generalized 
ChPT ?). 

5.4. m U: d.s to order a 3 from sum rules + 
ChPT 

One should note from Table 1 the consistency of 
the results from the pion channel and the one 
from the direct extraction of (ip^j which is an a 
posteriori support of the validity of the OPE for 
the 7r-sum rule in the working region, and signals 
the absence of the large effects due to instantons, 
which may break the OPE. Using the ratios from 
ChPT in section (||), one can deduce in units of 
MeV, the value of the running masses at 2 GeV 
to order a 3 given in Table 6. We have used the 
conversion factor: 

mi(l) ~ (1.38 ±0.06)771,(2), (20) 

for running, to order a 3 , the results from 1 to 
2 GeV, which corresponds to the average value 
of the QCD scale A 3 ~ (375 ± 50) MeV from 
PDG and others [Q . I remind that the errors 
in these determinations already take into account 
the systematics of the method (see Table 1). 

6. DIRECT EXTRACTIONS OF 777, 



Table 4 

Direct extractions of m s (l) in MeV to order a^. 
Channels m s (1) Sources Authors 



Kaon SR 165 ± 15 
155 ±25 
155 ± 25 

Scalar SR 203 ± 20 
143 ± 17 
160 ±30 
159 ± 11 
175 ± 48 

(r-like </> SR 173 ± 33 
e+-e _ data 176 ± 31 
±r-decay ) 186 ± 31 
179 ± 39 

AS = -1 2tet% 
part of 200 ± 50 
r-decay 164 ± 33 
213 ± 82 




Aio 
Ai 



SN89* ga 
DPS99 || 
Largest Range 

CPS97 gl[ 
CFNP971621 
J98 p 
M99 |64 
Largest Range 

SN95,99 



Largest Range 

Aleph99 
CKP98 
PP99 (70) 
Largest Range 



Average of Largest Ranges 166.7 ± 18.8 



6.1. Pseudoscalar sum rules 

In the strange quark channel, we quote in Table 
4 the results from and J59], and consider the 
largest range spanned by these previously quoted 
results. We consider that this conservative range 
already takes into account and may even overesti- 
mate the systematics of the method. One should 
notice here that, unlike the case of the pion, 
the result is less sensitive to the contribution of 
the higher states continuum due to the relatively 
higher value of M^, though the parametrization 
of the spectral function still gives larger errors 
than the QCD series. 




m s has 
4l by us- 



6.2. Scalar sum rules 

Following the pioneer's analysis of fl5 
been obtained by different authors 
ing the Kn phase shift data for parametrizing the 
spectral function. The different values obtained 
from this channel to order is given in Table 
4. Like in the case of the pseudoscalar channel, 
the errors are dominated by the uncertainties for 
parametrizing the spectral function. 

6.3. m s from e + e + r-decay data 

One can combine the e + e~ — * I = 0, 1 hadrons 
and the rotated recent AS = component of the 
r-decay data in order to extract m s . Unlike pre- 
vious sum rules, one has the advantage to have a 
complete measurement of the spectral function in 
the region covered by the analysis. We shall work 
with: 



Rt,4> = 



31 Kd 



-S. 



EW 



A./: 



ds 1 



M T 



1 



(21) 



2TTQ 2 

2s 

M?J M r 2 ' 

and the SV^-breaking combinations |37l , |05| : 

Ai^ = R T: i — R T ,<f>-, Aio = R T ,i — 3i? Tj o , (22) 

which vanish in the SU(3) symmetry limit; Aio 
involves the difference of the isoscalar (R T ,o) 
and isovector (R T ,i) sum rules a la Das-mathur- 
Okubo |36|. The PT series converges quite well 
at the optimization scale of about 1.6 GeV J65| . 
E.g, normalized to rn 2 sl one has: 



30.4af + (173.4 ± 109.2)4 (23) 



It has been argued in [mj that Aio can be af- 
fected by large SU(2) breakings, but this claim 
has not been confirmed from the result based on 
the other sum rules not affected by these terms 
|6q] . The largest range of values from different 
form of the sum rules is given in Table 4, which 
one can compare with the average of (178 ± 33) 
MeV given in |35| . An upper bound deduced from 
the positivity of i? r .0 is given in Table 3. 

6.4. m s from the AS = — 1 part of r-decay 

One can also extract m s from the Cabibbo sup- 
pressed channel of r-decay Jl4| , |68| -[70]| , using dif- 
ferent r-like moments. Unlike the case of the 



neutral 0-meson current, where the QCD series 
is more convergent, here the convergence is quite 
bad, such that one needs to select an appropri- 
ate combination (spin 1+0 pieces) for obtaining 
an acceptable result. Though a complete agree- 
ment has been obtained in the previous analysis 
of |6£j with the two other determinations [Qg8[ , 
a recent analysis in |70[] is lower and more pre- 
cise than the former, though still in agreement 
with the previous ones due to the generous er- 
rors given there. Ref. |7(]] argues that one should 
consider the previous results as an upper bound 
rather than a determination, an argument which 
needs to be confirmed. By inspecting the results 
in J7(| , we notice that the estimate decreases with 
increasing power of moments, rather than stabi- 
lizing. Therefore, it can be more appropriate to 
consider the conservative range spanned by the 
results from the three moments which is (180±68) 
MeV, rather than taking their average quoted in 
Table 4 from [^o|. This conservative value is in 
better agreement with the two other determina- 
tions. It is also interesting to notice that the re- 
sults from r-decay are in good agreement with 
the one from e + e~ data, an agreement which is 
a priori expected because of the similarity of the 
two approaches. 

6.5. m s and ra s j (m u + md) from QSSR 

More generally, one can notice from Table 4 that 
there is a total agreement of the sum rule results 
from different channels. As already mentioned, 
we expect that the largest ranges given in Table 
4 already include the QSSR systematics. Using 
the average of these largest ranges, one can de- 
duce the pure QSSR determination of m s , to or- 
der a' given in Table 4 and evolved to 2 GeV 
in Table 6. Combining this result with the sum 
rules determination of (m u + rrid) in Table 1 (sec- 
tion (5)), one can deduce pure QSSR prediction 
of m s /(m u + rrid) in Table 6, which is in agree- 
ment with and an independent test of the ChPT 
result 24.4 ± 1.5 given in section (3), though less 
accurate. 



7. QSSR + ChPT FINAL RESULTS 

We take the average of m s from m u + (Table 
1) + the ChPT ratio and from the direct deter- 
mination in Table 4. Then, we obtain the final 
average from QSSR+ChPT to order a, in Ta- 
ble 6. Combined with the ratios from ChPT in 
section 2, this value leads to the values of m Ui d 
given in Table 6. As already discussed in previ- 
ous sections, the quoted error already include the 
systematics of the methods. The size of the error 
is within the expected accuracy of the sum rule 
results. Using Eq. 2, it is trivial to extract the 
value of the invariant mass m;. One obtains in 
units of MeV: 

m u = 3.9 ±0.7, rh d = 7.1 ± 0.8 

m s = 133.3 ±18.8, (24) 

where the error is larger than the corresponding 
running mass due to A in the evolution procedure. 

8. COMPARISON WITH THE LATTICE 

8.1. Lattice approaches for/by non-experts 

One usually starts from the QCD action and par- 
tition function: 

Z = J VA^ det M J d ^(-i^G^) ( 25 ) 

integrated over gauge field configurations. The 
fermion contributions are included into the non- 
local det M term. For the analysis, one works like 
in the sum rule approach, with the 2-point cor- 
relator defined in previous sections, which is sat- 
urated by the intermediate states |n)(n|. In this 
way, the two-point correlator can be expressed as: 

^<0|J(o:)(J(0)) t |0) = 

X 

^{0\J( x )\n){n\(J(0))U0)^ (26) 

where the zero momentum states E n tend to the 
masses M n of the resonances. In the (ideal) 
asymptotic limit t — > oo, the exponential fac- 
tor kills the effect of the different excitations, 
such that the lowest ground state contribution 
dominates. In practice, this approximation is ex- 
pected to be realized when the splitting between 



the ground state and the radial excitation is large 
enough. 

8.2. Practical limits of the lattice 

Besides the usual statistical and hnite size 
(about 1% if the lattice size L > 3 fermi and 
m-nL > 6), errors inherent to the lattice, which 
can be minimized using modern technology, there 
are still large uncertainties related to the uses of 
field theory on the lattice due to the finite values 
of the lattice spacing a: 

• The different operators mix at finite a. 

• The discretization errors specific to each ac- 
tions, which are 0(a) for the Wilson (explicit 
breaking of chiral symmetry (x^)) and Domain 
wall (extra 5th dimension in order to preserve 
XS) actions, 0(a 2 ) for the staggered (reduction 
of quark couplings with high-momenta gluons) 
and 0(aa s ) for the Clover (inclusion of the mixed 
quark-gluon operator) actions. For typical values 
of 1/a ~ 2 GeV, the error is « 10-30%, which can 
be reduced by computing at different values of a. 

• The well-known quenched approximation (no 
inclusion of the fermion contribution In Det M), 
which implies a modification of x>S with unphys- 
ical singularities for m q — or practically for 
m q < 77i s /3 (recall that in this approximation: 
M v i w inn = (= large A c -limit)), which in- 
duces an error of about 20% that can be esti- 
mated from the deviation of the predictions from 
the observed meson masses and couplings or / and 
from the choices of the mesons for setting the 
scale (string tension). 

• The extrapolation of the results to light quark 
masses with the help of the meson mass depen- 
dence expected from ChPT, which for typical val- 
ues 1/a — 2 GeV, and keeping m T L > 6, one re- 
quires L/a > 90 in order to avoid finite volume 
effects. At present, L/a ss 32 (quenched) and 
L/a w 24 (unquenched) far below this limit. 

• The errors due to the matching of the lattice 
and the continuum at a typical lattice conversion 
scale of 2 GeV can be minimized using the non- 
perturbative renormalization. 

8.3. Lattice results and estimated errors 

From the previous discussions, we consider that: 

• The conservative quenched lattice errors are 



Table 5 

Lattice Quenched Results up to NNLO since 98. 



Group 


m s (2) 


Sources 


Comments 


98 








OHIojf 


J 129± 23 




Staggered 


ape frc 


J 130±18 


K* 


NLO, AWI 








Wilson, Clover 




121±13 


K, <t> 


NLO 








NPR+AWI 


99 








CP-PACS 


H 143±6 


4> 


AWI+VWI 








Wilson 




115±2 


K 


qChPT 


jlqcdJ; 


'J 129±12 


4> 


AWI±VWI 








Kogut-Suss. 




106±7 


K 


qChPT 


a-UKQCD[F| 97±4 


!k 


NNLO 






N 


/ 10% 


DESYfra 105±4 


K* 


NNLO 








AWI 








I. Wilson, 


BNLjf? 


] 95±26 


U-K 


Domain 


BNLjT? 


130 ±21 




Walls 


APE[£E 


| 114±9 


Q-Prop. 


NNLO 


Average 112.9 


±1.5 (stat.; 


±22.3 (syst.) 



about 20%. 

• The extraction of m Uy d is less reliable than m s . 
Therefore, we shall only consider the value of m s 
obtained from the lattice which we shall compare 
with the one obtained in previous sections. Lat- 
tice results prior 98 have been already reviewed in 
fn|. The different results for 98 and 99 are given 
in Table 5 for different actions, where one can see 
a large spread of predictions, which with the given 
errors are inconsistent each others. We mainly at- 
tribute the source of this discrepancy to the un- 
derestimate of the systematic errors given there. 
Most of these results have been obtained using 
the non-perturbative renormalization fl§0| , and 
Ward identities for the axial (AWI) and/or vec- 



Table 6 

Summary for m„ )( j iS (2) in MeV 



Sources m s 2m s /(m u + md) m u + mj, m u /m,d m u rrid 



QSSR 

Table 1 
Table 4 
Average 

LATTICE 

Table 5: quenched 
nf — 2 dynamical Q 



115.8 ± 19.7 

120.8 ± 14.6 

118.9 ± 12.2 



112.9 ± 22.4 * 

97 ill-v/2 * 



GRAND AVERAGE 110.9 ± 8.8 

QSSR BOUNDS 90 < ... < 168 

Tables 2 and 3 



$= ChPT+ 
25.5 ± 4.8 

f ChPT =*> 



-ChPT 



9.5 ± 1.4 +ChPT 
9.8 ± 1.2 +ChPT 

less reliable 
less reliable 

9.1 ± 0.9 +ChPT 

7 < ... < 12 



3.4 ±0.6 6.1 ±0.8 

3.5 ±0.4 6.3 ±0.8 



3.2 ±0.3 5.9 ±0.5 



* Our conservative error estimate 



tor (VWI) currents, and constraints from ChPT 
in the extrapolation procedure. There is also a 
systematic discrepancy between the results from 
the kaon and <f>, K* or the nucleon masses, which 
might be some indications of the quenching er- 
rors, while e.g. the splitting of the K-K* is no 
longer resolved. As one can see from Table 5, the 
quenched lattice predictions are in the range: 

m s (2) ~ (69 — 181) MeV , (27) 

where part of this range is already excluded by the 
bounds given in Tables 2 and 3. Instead, one can 
also quote (to be taken carefully) the naive aver- 
age given in Table 5 at NNLO Q where we have 
added our guessed 20% estimate of the lattice sys- 
tematic errors based on the previous comments. 
At this approximation, where a comparison with 
the previous results from QCD spectral sum rules 
is meaningful, one can notice a surprisingly good 
agreement. 

Some attempts to put dynamical fermions have 

7 However, as discussed in |79| , the inclusion of the higher 
order a s corrections obtained in M] in the conversion 
of the lattice to the continuum results tends to decrease 
slightly the value of m s by about 3%. 



been done in ]82[| , and more recently with 2 
flavours in Q. In |83], some of the problems 
encountered in the quenched approximation (dis- 
crepancy between the </> and K results,...) seems 
to be resolved. Though promising, the approach 
is not enough mature for the different systemat- 
ics to be fully under control. We quote in Table 
6 this result adopting a more conservative error 
than the original one. 

9. SUMMARY 

We have reviewed the different determinations 
of light quark masses from ChPT, QCD spectral 
sum rules (QSSR), e + e~ and T-decay data, and 
compared the one of the strange quark mass with 
the recent lattice results: 

• The sum of light quark masses m u + to order 

from different QSSR analysis is given in Table 
1 and the resulting average value. 

• Lower (resp. upper) bounds based on the pos- 
itivity and analyticity properties of the spectral 
functions are given in Table 2 (resp. Table 3). 

• Different direct determinations of the strange 
quark mass to order are compiled in Table 4. 

• Combined results from these four methods lead 



to the final average given in Table 6 to order af, 
where the errors are typically the 10% systemat- 
ics of the QSSR approach. An eventual failure 
of this result should signal a new phenomena not 
accounted for in the OPE discussed in this paper. 

• We have compared this final result with the re- 
cent (after 98) lattice determinations (Table 5) 
which belong in the range given by Eq. 27 and 
which lead to the average in Table 5. 

• Within the present uncertainties of various ap- 
proaches, we consider that there is a good agree- 
ment between the previous sum rule and lattice 
results. Attempting to give the final QCD value, 
we average the different results in Table 6, and 
deduce the QCD Grand Average given in this ta- 
ble, to be used carefully. 

• However, we expect that a future high-precision 
measurement of the light quark masses will be 
difficult to reach due to the systematic errors in- 
herent to each methods, which, often, different 
authors do not include into their results ! 

• Finally, one should remind that, in the phe- 
nomenological analysis, one should use the value 
of m s into the expression of any hadronic matrix 
elements or / and observables which are known at 
the same level of approximation. This consistency 
condition is not often respected in the literature. 

10. APPLICATION TO e'/e 

One of the most fashionable applications of the 
previous result is the one to the CP violating 
parameters e'/e, where e 1 is related to A\Kj, — > 
(tttt) 1-2] /A[Ks — > (7r7r)j=o] and characterizes the 
(direct) CP-violation in the decay amplitude of 
K -> tttt; e = A[K L -> {mr) I=Q ]/A[K s -> 
(7T7r)/=o] is the (indirect) CP-violation from K°- 
K° mixing. It is known || that the dominant 
effects in the analysis of Ki,,s ~ * ( 7r7r )/=o,2 am- 
plitudes are due to the QCD and electroweak pen- 
guin operators: 



r>3/2 , 2 

B s /m s 



0(1/N) 



(28) 



Q 6 = (s a d ) v _ A 2J tyartp) 



V+A 



Bl /2 /m 2 s +0(l/N) 

o (s a d p ) v _ A (^i> a ip(3) v+A 



For M t ^ 165 GeV and A 4 ~ 340 MeV, the sim- 
plified SM prediction without (resp. with) the in- 
clusion of final-state interaction effects 84 1, is||: 



9.75(resp. 15.34)ImA t 



110 MeV 

m s (2) 



B, 



1/2 



0.54(resp.0.32)Bg /2 



(29) 



where At — VtdV t * is expressed in terms of the 
CKM matrix elements. Using Im At ~ (1.34 ± 
0.30) x 10~ 4 IH, and the measured value §] 
W/e)ex P — (21.4 ± 4.0) x 10~ 4 , one can deduce, 
from the average value and the lower bound of 
m s in Table 6 f\ 



bI' 2 - 0.54(0.32)5, 



3/2 
8 



~ 1.6 ± 0.4(1.1 ± 0.3) , 
> 1.1 ±0.2(0.7 ±0.1), 
(30) 



which signals a violation of about 1 *~ 3er for 
the leading 1/N vacuum saturation prediction ^] : 
B^ 2 w B^ 2 sa 1. This result and the final-state 
interaction effects should be tested using more ac- 
curate non-perturbative calculations. It is only 
after performing these tests that one can make a 
sharper conclusion on the SM prediction of the 
CP-violation 0. 
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QUESTIONS AND DISCUSSIONS 

Lively general discussions, many questions and 
some comments have followed this talk, As they 
have been also addressed to the previous talks in 
this session, they have not been reported here. 



